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Abstract — In the literature the performance of quantum data 
transmission systems is usually evaluated in the absence of 
thermal noise. A more realistic approach taking into account the 
thermal noise is intrinsically more difficult because it requires 
dealing with Glauber coherent states in an infinite-dimensional 
space. In particular, the exact evaluation of the optimal mea- 
surement operators is a very difficult task, and numerical 
approximation is unavoidable. The paper faces the problem by 
approximating the P-representation of the noisy quantum states 
with a large but finite numbers of terms and applying to them 
the square root measurement (SRM) approach. Comparisons 
with cases where the exact solution are known show that the 
SRM approach gives quite accurate results. As application, the 
performance of quadrature amplitude modulation (QAM) and 
phase shift keying (PSK) systems is considered. In spite of the 
fact that the SRM approach is not optimal and overestimates 
the error probability, also in these cases the quantum detection 
maintains its superiority with respect to the classical homodyne 
detection. 

Index Terms — Quantum detection, square root measurement, 
geometrically uniform states, thermal noise, quadrature ampli- 
tude modulation (QAM), phase shift keying (PSK). 



I. Introduction 

Transmission of information through a quantum channel is 
mainly affected by an uncertainty which is intrinsically related 
to the quantum mechanics laws. In the language of classical 
optical systems, this uncertainty corresponds to the so called 
shot noise. Another reason of uncertainty is the presence 
of tliennal noise, as in classical systems. Beginning from 
1970's, a lot of research work has been devoted to the 
quantum detection problem [10], that may be summarized in 
the following terms. The transmitter sends a quantum signal 
through a quantum channel, which forces the receiver to 
assume one among a finite number of states. The detector tries 
to guess the state by an adequate set of quantum measurements 
and the problem arises of finding the measurement set which 
optimizes the detection, according to some predefined fidelity 
criterion (usually the minimum error probability). Necessary 
and sufficient conditions for the optimal measurement set have 
been found in pioneering papers by Holevo [11] and by Yuen 
et at. [15]. 

Unfortunately, even though the optimal measurement set is 
completely characterized, analytical closed-form solutions for 
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the measurement set, indeed for the minimum error prob- 
ability, are not available, in general. Then, it is needed to 
resort to a numerical evaluation based on convex semidefinite 
programming [6]. However, under some rotation symmetry 
constraint, a simple measurement, introduced by Hausladen et 
at. [8] and known as square root measurement (SRM), turns 
out to be optimum. The SRM has the remarkable advantage 
that it is straightforwardly evaluated starting from the possible 
states. Moreover, also when it is not optimal, the SRM often 
gives "pretty good" upper bounds on the performance of 
optimal detectors. 

Our paper starts just from these important results on SRM 
for studying quantum data transmission systems in the pres- 
ence of thermal noise. After the pioneering work in [10], the 
problem of quantum detection in a noisy environment has 
received scarce attention in the literature. To the best of our 
knowledge, only a correspondence by Sasaki et al. [13] on 
quantum on-off keying (OOK) and a technical report [14] 
on quantum PSK attempt to afford an approximate analysis 
of quantum detection of coherent states. This delay is due 
to the difficulties of performing efficient approximations in 
the numerical performance evaluation. On the other hand, 
exploiting the new perspectives open by the extension of the 
SRM to mixed states [6], we apply this approach to a quantum 
noisy channel according to the Glauber theory on coherent 
states. 

The paper is organized as follows. In Section II we review 
quantum detection fundamentals. In Sections III the SRM tech- 
niques are recalled and in Section IV the key problem of the 
finite-dimensional factorization of the Glauber representation 
of noisy states is discussed. Finally, in Sections V and VI the 
SRM approach is applied to QAM and PSK quantum systems 
(the same systems considered by Kato et al. [12] in the absence 
of thermal noise). 

II. Quantum Data Transmission System 

In this section we recall some basic facts about quantum 
detection following the scheme of Fig. 1. For a detailed 
treatment the reader is referred to [10] and for a more recent 
survey to [2]. 

pure state density op. 
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Fig. 1 — Model of a quantum data transmission system. 
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A. General Model 

A classical source emits a symbol a drawn from a finite 
alphabet yi = {0,l,...,m— 1} with a given prior probability 
distribution qi = V[a = i] , i = 0, . . . , m — 1. On the basis 
of the symbol a emitted by the source, the transmitter sends 
a quantum state \^a) through a quantum channel (e.g., a laser 
pulse through an optical fiber). As a consequence, the received 
state is one of m possible quantum states and the detection 
device performs a set of measurements in order to guess the 
received state and consequently the original symbol. 

In the ideal case, i.e. neglecting thermal noise, a set of 
TO pure states is seen by the receiver, which are a replica 
of the transmitted states |7a). In the presence of thermal 
noise the received states become noisy (or mixed) and are 
described by a set of density operators pi, f = 0, . . . , to — 1, 
which are Hermitian, positive semidefinite (PSD) and have unit 
trace, Tr(pj) = 1. The description through density operators 
represents the general case, since it is comprehensive of the 
pure state case, in which pi reduces to the rank-one operator 
Pi = |7i)(7i|- 



B. Quantum Detection Theory. Available Results. 

Quantum theory postulates that a detection device for 
choosing among the possible states is given by a positive 
operator valued measurement (POVM), i.e. a set of m oper- 
ators Ho, ... , Ilm-i that are Hermitian, PSD and resolve the 
identity operator of the Hilbert space J{, namely X^I^o^ — 
I-K- Then, the probabiUty that the detection system reveals the 
state j, provided that the state density operator is pi, is given 
by 

p{j\i) =Tr:{pin.j) , i,j = 0,...,m-l. (1) 
In particular, the probability of correct detection becomes 

m— 1 m— 1 

Pc=J2qi P{i\i) = ^ % Tr{piUi) . (2) 

i=0 i=0 

For pure states, that is with pi = rank-one POVMs 

of the form Uj = |Atj)(/Zj| can be used, where \iJ.j) are called 
measurement vectors. Then, (1) reduces to = |(7i|Atj)p. 

The optimization of the detection scheme reduces to finding 
the POVM operators 11^ that maximize J = X)™ TT{qiPiIli) 
under the constraints that the Ilj are PSD and resolve the 
identity I^. The maximum of J is the optimal probability 
of correct detection. Clearly, this is a convex semidefinite 
programming problem in the real space of the Hermitian 
operators. 

On the other hand, analytical closed-form results are avail- 
able only for the particular class of pure states exhibiting 
the so called geometrically uniform symmetry (GUS). In this 
case Ban et al. [1] have shown that the optimal POVMs 
are given by the SRM. This particular solution has been 
thoroughly discussed by Eldar and Forney [4]. Recently the 
SRM approach has been extended to mixed states by Eldar et 
al. [6]. 



III. The SRM Techniques 

In this section we review the SRM techniques, having in 
mind that our final application will be the optical quantum 
transmissions, where quantum states and density operators 
should be formulated according to the Glauber theory (see 
Section IV). The SRM technique is here considered in the 
general case of mixed states, following Eldar, Megretski and 
Verghese [6]. 

A. General Formulation 

We start from a constellation of m density operators 
po, . . . , pm-i in an n-dimensional Hilbert space Jf. The key 
of the SRM approach is the factorization of each density 
operator in the form pi = 7^ 7* for some complex matrices 
7i, e.g., via the eigendecomposition of pi. The factorization 
is not unique, but the ambiguity is irrelevant for the quantum 
decision. If pi has rank < n, the factor 7, can be chosen to 
have dimensions n x r^. In [6] 7^ is referred as a factor of pi, 
but, more specifically, we call 7j a state factor. Since the i-th 
optimal measurement operator can be chosen with rank not 
greater than the rank of pi [5], the search can be confined to 
POVMs of the form 11^ = pi p* where pi are nx ri complex 
matrices. We refer to pi as measurement factors. 

The (generalized) state and the measurement matrices are 
obtained by storing the corresponding n x rj factors as block- 
colunrn vectors, namely F = [70,71, ... ,7m-i] and M = 
[/io, /ii, • • ■ , Pm-i] - The dimensions of both F and M aienxk 
with k = ro-\ h Vm-i- 

From the state matrix F we form the k x k Gram matrix 
G = F* F and also the n x n matrix T = FF*, sometimes 
called Gram operator. 

At this point the SRM method is used to provide the mea- 
surement matrix M. The first step is the eigendecompositions 
of G and T, namely G = VAqV*, T = UAtU*, where 
U, V are unitary and Aq, At are diagonal. Note that G and 
T are both PSD with the same rank r and have the same 
set of positive eigenvalues [4]. From the eigendecompositions 
we can find the inverse of the square root of G and T as 
= VA'^^^V* and T-^/^ = UA'^'^U*, where 
the inverses must be intended in the generalized Moore- 
Penrose sense [4]. Finally, the measurement matrix is given 
by M = T^i/^p. 

An alternative evaluation of M is obtained through the 
Gram matrix. In fact, it can be shown [4] using the singular- 
value decomposition, that the measurement matrix is also 
given by M = r G~^/^. Then, the evaluation of the transition 
probabilities follows from (1), namely 

P{j\i) =Tr(Pinj) = Tr(7i7>j^*) 

=Tr(/x*7i7>i)=Tr(%S;j ^> 

where Bji is the (j, z)-th block of the matrix M*T = 
(0-i/2p*p ^ (^1/2 rpj^gjj^ foj. jj^g evaluation of p{j\i) we 

have to partition G^/^ into blocks. Finally, the probabiUty of 
correct decision becomes 

^ m— 1 
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This twofold possibility, via T or G^^/^, is very im- 
portant for an efficient computation, particularly when k < n. 

B. SRM with Geometrically Uniform Symmetry 

The SRM is simplified and provides pecuUar properties 

if the state constellation exhibits the geometrically uniform 
symmetry (GUS), that is if there exists a unitary operator' S, 
such that 5™ = 7j{ and pi = S^poS~\ The operator S and 
the density operator po are said the generating operator and 
the generating density of the constellation, respectively. For 
the mixed states the factorization po = 7o7o leads to the form 
7i = S"7o- Note that with GUS all the state factors 7, have 
the same rank h as 70. 

Since 5 is unitary, its eigendecomposition has the form S = 
YAY* = Yljlg ^j\yj){yj\, where F is a unitary matrix of 
order n and the eigenvalues Xj collected in the diagonal matrix 
A have unit amplitude. Moreover, because of = /jc, the 
diagonal matrix has the form A = diag [VF,^", . . . , Wm ^^] 
where Wm = e*^'^/'" and the exponents ri are suitable integers 
with < rj < m. By collecting the terms with equal eigenval- 
ues in the eigendecomposition, one gets S = Y^'^Sq W^Y^, 
where are projector operators, so that YhY^ = Y^dhk- 

The Gram matrix G = r*r, of order mh, is formed by the 
blocks of order h 



m— 1 



m— 1 



where = m^Q Y^ 70. Since Grs depends only on the 
difference r — s mod m the Gram matrix is block circulant. 
But, the important point is that this property yields an explicit 
decomposition for G, namely 



G = V^.hDV* 



(4) 



where D = diag {-Do, • • • , -Dm-i} and Vm,h is the hm x hm 
matrix Vm,h = IKVV^)^™ -^'Jl with Ih identity matrix 
of order h. As a consequence, the diagonal blocks of the 
matrix are given by the Discrete Fourier Transform (DFT) 
of the first "block row" of the Gram matrix G, namely 

Note that (4) is not a standard eigendecomposition because 
the blocks Dj are not diagonal matrices. To find the square root 
of G we have to evaluate the square root of D with an eigende- 
composition of each block Dj. Since these are PSD Hermitian 
square matrices, their square roots matrices D^^^'^ can be 
calculated to construct D^^^'^ = di&g [D^^^^^ , . . . , D^^/^]. 
Finally, we obtain G±V2 = Vm,h ^^^1"^ V^^u whose (r, s) 
block is given by 

^ m— 1 

[^^'^'1- = ^ Y.'^^^-'^'^T" ■ (5) 

V '''' n 

Now, the probabilities are obtained by applying (3) 

with Bji = (G^l'^)ji. In particular, ■p(i\i) are independent of 

'The GUS can be generalized over a group of unitary operators [6], but in 
our applications this generalization is not needed. 



i and give the probabiUty of correct detection. The expUcit 
result is 



:p(z|i)=Tr[i?2] =Tr[(Gi/2) 

- m — 1 



fe=0 



(6) 



Finally, we recall that the SRM is optimal for GUS pure 
states, but not for GUS mixed states, at least in general. In 
particular, the sufficient condition for optimaUty given in [6] 
fails in all our examples of appUcation. 

IV. Signal and Noise in Quantum Optical 
Communications 

In this section we recall the quantum environment for the sig- 
nal and thermal noise in optical communications. The correct 
settlement is provided by the celebrated Glauber theory [7], 
which represents signal and noise in an infinite dimensional 
Hilbert space. 

A. Representation of Coherent States 

The quantum model of a coherent state representing a 
monochromatic electromagnetic radiation produced by a laser 
is formulated in an infinite dimensional Hilbert space !K 
equipped with an orthonormal basis |n), n = 0, 1, 2, . . ., where 
|n) are called number eigenstates. Each state \n) is said 
to contain exactly n photons. In this context the Glauber 
representation of a single radiation mode is given by the ket 



I7) = e 



(7) 



where 7 is the complex envelope that specifies the mode. Thus, 
for each 7 e C a coherent state (or Glauber state) is defined; 
in particular, the state |0) obtained with 7 = represents the 
ground state. The probabiUty of obtaining exactly m photons 
is governed by the Poisson distribution p{m\^) = |(to|7)P = 
exp(-|7p) |7p'"/m! with mean |7p. Hence, = |7p 
represents the average number of photons when the system 
is in the coherent state I7). We recall that the Glauber states 
are not orthogonal, since the inner product of two Glauber 
states is given by 



-i(|a|^ + |/3|^-2a*/3) 



(8) 



The representation (7) is only valid when the receiver 
observes a pure state with a known parameter 7, which in 
the context of communications may be regarded as the signal. 
In the presence of thermal (or background) noise the signal 7 
becomes uncertain and must be represented through a density 
operator The Glauber theory [7][10] states that the density 
operator is given by 



(9) 



that is by a continuous mixtore of coherent states. In (9) 
the parameter K represents the average number of photons 
associated with the thermal noise; it is given by K = 
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l/[exp (hiy/kTo) — !], with h Planck's constant, k Boltzmann's 
constant, i/ optical frequency and Tq absolute temperature of 
the receiver. Hence, the representation of "signal plus noise" 
depends only on the two parameters: 1) 7 € C, which deter- 
mines the nominal coherent state I7), and 2) 3\f representing 
the average number of noise photons. When Dsf = 0, that is 
in the absence of noise, relation (9) degenerates into the pure 
state density operator ^(7) = |7)(7|. 

B. Discretization of the Density Operators 

An infinite matrix representation ||pmn|| of the density 
operator (9) is obtained in terms of the orthonormal basis of 
the number eigenstates \n), namely pmn = {'ni\p{l)\n) , and 
the expression of the m n entry is [9] 



Pmn (7) =(1 




-(1-,;)|7| = 



jn- 



- J^) and L: 



(10) 



where 0<m<n, jy^O, v = J^/{1 + J^) and L;^-'"(x) are 
the generalized Laguerre polynomials. The entries for m > n 
are obtained by the symmetry Pnmil) = Pmnil)- The matrix 
is infinite dimensional and not diagonal. The diagonal elements 
Pmrnil) give the probabilities of obtaining exactly m photons 
when the quantum system is in the noisy state ^(7) [9]. From 
(10) we have 



PL{m) =Pmm{l) 

=(1 - v)?;™ e-^i-")^^ L„((l - vfN^/v 



(11) 



which represents the Laguerre distribution (L„i{x) = i°„(x) 
are the ordinary Laguerre polynomials). The mean and the 
variance of distribution (11) are and + 2N^'N + 

^^(Jvf + 1), respectively. For the ground state I7) = |0) 
the above expressions degenerate. The matrix representation 
becomes diagonal, namely 



Pmn{0) = <5mn(l - v)v'^ 



(12) 



and 



the distribution becomes geometrical: pcijn) = 
Pmm{0) = (1 ~ v)v"^. The infinite dimension matrix 
||Pmn(7)||)0 < m,n < 00 gives a correct representation of 
the density operator. But, for the SRM, which is based on 



eigendecompositions, we need a finite dimensional approxima- 
tion by a truncation to iVe terms. For the choice of to get 
a given accuracy we follow the quasi-unitary trace criterion, 
which is based on the fact that a density operator has a unitary 
trace. Then, we choose as the smallest integer such that 



E 

m=0 



= E 

m—O 



(7) = Yl PL{m) > 1-e 



where e is the required accuracy. Thus, for a given e, A^^ can 
be evaluated using the Laguerre distribution (11). 

C. Factorization of the Density Operators 

Once established the finite A^^ x approximation of 
the density operators, for the SRM we need a factorization 
of the form p{a) = 7(Q!)7*(a) for a convenient matrix 
7(a), which we call state factor This is obtained from the 
eigendecomposition of p{a), namely 



i=l 



UrKu: 



where r is the rank of p{a), Ur is x r and collects the 
eigenvectors \ui) corresponding to the r positive eigenvalues 
Af , which are assumed in a decreasing order, and is r x r 
diagonal collecting the Af. Hence, 7(a) = UrAr is a correct 
factor of p(a). (For a = (ground state) the factorization 
is immediate since p{0) is diagonal and from (12) we find 

7(0) = Vp(0) = ||5mnV(l-yKII-) 

A critical point in the numerical evaluation is the choice of 
the rank r, given by the number of the numerically relevant 
positive eigenvalues. To clarify the problem we develop a 
specific case: a = \/5 , iVc = 5 , K = 0.1 . f = 10^^ —^ 
= 20. Now, in theory p{a) has a full rank r = N^, as we 
can see from the list of its eigenvalues obtained with a great 
accuracy 



0.150285 , 0.00231095 , 0.0000353779 , 5.2072510" 

-19 



, 7.2487410" 
-4.1356410"" 



9.4515710"" , 1.1460310"^^ , . . . , 1.344510" 
- 3.1086710"^^ , 2.3218610"^^ , -1.3763110"^^ , 3.0177510"^^ 

but in practice we can limit to take only the first 3 eigenvalues, 
neglecting the remaining, which means that we assume as a 
"practical" rank r = 3. As a check, the reconstruction of p{a) 
obtained in such a way assures an accuracy < 10"^. To find 
the "practical" rank in general we consider the reconstruction 



Table 1 ; Values of and Ni, of the dimensions of state factors for some values of the average 
number of photons N~f and of the thermal noise parameter X 
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error Ap = p — 7r7*> where the factor 7^ is obtained by 
considering only r eigenvalues. Then, we can evaluate the 
maximum error or the mean square error (m.s.e.) as a function 
of r and we choose r = to achieve a given accuracy v. 

In Table 1 we give a collection of Af^ and Ni, for some values 
of and K, obtained with the accuracies e = v = 10~^. 
was obtained considering the m.s.e. 



V. Application to QAM Modulation 

The m-QAM constellation is defined starting from the auxil- 
iary alphabet = {-{L - 1) + 2(i - 1)| i = 1, 2, . . . , i} 
with L = 2, 3, . . . and is given by the m Glauber states 



|7„„) = |A(u + iu)) , u,veAL 



(13) 



with m = L"^. This constellation has not the GUS and 
therefore the SRM must be applied in the general form. In 
(13) A is a scale factor, which determines the average number 
of signal photons, specifically 



1)A2 



2, 



1)A^ 



(14) 



For instance, for the 16-QAM we find = lOA^. 

In the case of pure states the first step is the evaluation of 
the Gram matrix G, whose elements are given by the inner 
products^ 

{luv\lu'v') = (A(u + iz;)|A(u' + w')) 

= exp{-iA2[(M' - uf + {v' - vf - 2i{u'v - v'u)]} 
u, V, u', v' e Al ■ 

Then, the eigendecomposition G = VAgV*, the evaluation 
of G^/^ and of the probabilities can be carried out without 
approximation and with a low computational complexity since 
the dimensions involved are only m x m [12]. 

A. Application of SRM in the Presence of Noise 

In the presence of thermal noise the only problem is the 
management of approximations since the density operators 
must be approximated by finite-dimensional matrices, as dis- 
cussed in Section IV. In the QAM format the average number 
of photons = |(m + w)Ap, is not uniform over the 
constellation, varying from = 2A^ for the inner symbols 
to = 2{L — 1)^A^ for the corner symbols. The reduced 
dimensions of the Hilbert space n = must be chosen 
considering the maximum N^^max. = 2{L — 1)^A^. Then, 
assuming Ng as fundamental parameter, for the choice of N^: 
we have to consider that 



N. 



7max 



= 3[(L - 1) V(L2 _ 1)] Ns = 3[(L - 1)/{L + 1)] iV, 



For instance, for 16-QAM we find A^T-max = l-SNg. 

We sketch an example to illustrate the dimensions involved 
in the 16-QAM. With iV^ = 4 and ?^ = 0.1 we find A^^- 



7max 



7.2 and we choose = 40 assuring an accuracy e = 10 '. 
The dimensions of the pi are 40 x 40 and they are factored into 



matrices 7^ of dimensions 40 x 8. The dimensions of F, T and 
G are 40 X 128, 40 x 40, 128 x 128, respectively. So, it is more 
efficient to compute T~^/^ rather than G='=^/^. With the above 
choices we find the following diagonal transition probabilities: 
0.875749 for inner states, 0.916501 for side states, 0.947767 
for comer states, and the average error probability is Pg = 
0.08587. 

The SRM approach has been applied systematically to 
evaluate the error probabihty Pg in the 16-QAM systems 
following the steps outlined above. The results are illustrated 
in Fig. 2, where P,, is plotted versus the average number 
of photons per symbol Ns for some values of the thermal 
noise parameter DnT. In particular, the curve for K = 0, which 
refers to the absence of thermal noise, was checked with the 
results obtained with pure states and a perfect agreement has 
been found. To assure an overall accuracy of e = lO"'^, the 
dimensions of the Hilbert space have been set to the value 
n = AT, = 130. 

In Fig. 2 the performance of the SRM quantum detector 
is also compared with the performance of the classical ho- 
modyne detector, for which a Gaussian additive model with 
SRN=4A^s/(l + 21^) results.-' In the absence of thermal noise 
an improvement of about 3 dB over homodyne detection 
is confirmed. As it was expected, this improvement rapidly 
reduces as thermal noise increases. (For a comparison with 
optimum detection see the end of Section VI). 
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Fig. 2 — Error probability in 16-QAM versus Ns for some values 
of 3Nf. Solid lines refer to quantum detection ans dashed 
lines to classical homodyne detection. 



VI. Application to the PSK Modulation 

The constellation of a coherent PSK modulation format is 

given by the Glauber states |7i) = \ joW^),i = 0, . . . , m — 1, 
where, without loss of generahty, we assume that 70 is real 
positive. This is a GUS constellation with initial state I70) and 



depends on the four indexes u, v, u' , v' , but it can be arranged as an 
ordinary matrix using the lexicographic order for the indexes. 



^In Kato et al. [12] an analogous comparison is made with heterodyne 
detector. 
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generating operator 



The matrices Lfe are 



S = Y,W-\n){n\. 



(15) 



n=0 



With pure states the performance evaluation starts from the 
inner products Gqs = (70 lis), which can be obtained from 
(8), namely 

Go« =cxp[-72(l-W^4)], s = 0,l,...,m-l. (16) 

Then, the eigenvalues of the Gram matrix are obtained as the 
DFT of Gos 



m—1 



s=0 



Finally, the minimum error probability is 



^ m— 1 2 



(17) 



(18) 



fc=0 



The above expressions are obtained without numerical ap- 
proximations. Pe is a function of the alphabet size m and of 
the parameter 70, whose square A^^ = 70 gives the average 
number of photons per symbol. 

A. Application of the SRM-GUS in the Presence of Noise 

The T77 Glauber density operators obtained with a PSK 
constellation verify the GUS with generating operator S given 
by (15), which has infinite dimensions. In the SRM S must 
be reduced to a finite dimension n and then, in matrix form, 
it becomes 



5 = diag[iy^,,fc = 0,l,. 



. ,n- 



Its eigendecomposition S = Y AY* is trivial with F = /„ 
and A = S, and the matrices Lfe are given by 



Lk = m diag [Sk 



0,1, 



1] 



Now, for a given n, Ng ~ 79 and J^, the application of 
the SRM-GUS proceeds as follows: 1) evaluate the reference 
density operator po = Po(7o) from (10); 2) find the factor 
70 of po; 3) evaluate the blocks Dk = jnLkjo and find 

1 /2 

the square roots Z)^' by eigendecomposition; 4) evaluate the 
blocks Bii = {G^/'^)ii from (6a); 5) evaluate Pc from (6). 

We give a detailed example of calculation in the case of small 
dimensions (for reason of space). We consider the 4-PSK with 
Ng = \, = Q.\ and N^^ = 8, which assures an accuracy of 
e = 10~^. The reference density operator is the 8 x 8 matrix 



0.366 
0.333 
0.214 
0.112 
0.051 
0.021 
0.008 
0.003 



0.333 
0.336 
0.237 
0.136 
0.067 



0.214 
0.237 
0.183 
0.114 
0.060 



0.029 0.028 
0.012 0.012 

0.004 0.005 



0.112 
0.136 
0.114 
0.076 
0.044 
0.022 
0.010 
0.004 



0.051 
0.067 
0.060 
0.044 
0.027 
0.014 
0.007 
0.003 



0.021 
0.029 
0.028 
0.022 
0.014 
0.008 
0.004 
0.002 



0.008 
0.012 
0.012 
0.010 
0.007 
0.004 
0.002 
0.001 



0.003 
0.004 
0.005 
0.004 
0.003 
0.002 
0.001 
0.001 



Its practical rank is 5. From the eigendecomposition of po we 
obtain the 8x5 factor 



70 



-0.289 
-0.289 
-0.204 
-0.118 
-0.059 
-0.026 
-0.011 
-0.004 



0.087 
0.000 
-0.062 
-0.071 
-0.053 
-0.032 
-0.016 
-0.007 



-0.019 
0.019 
0.013 
-0.008 
-0.019 
-0.019 
-0.013 
-0.008 



-0.003 
0.006 
-0.002 
-0.005 
0.000 
0.004 
0.006 
0.005 



0.000 
0.001 
-0.002 
0.000 
0.001 
0.000 
-0.001 
-0.002 



40000000 
00000000 
00000000 









4 



00000000 
00000000 
00000000 



00000000 
04000000 
00000000 
00000000 
00000000 
00000400 
00000000 
00000000 



Then the evaluation of Dk and its square root, e.g. for = 0, 



gives 



Do 



r,l/2 _ 

^0 — 



0.348 -0.0S8 0.026 0.004 0.000 

-0.088 0.042 -0.002 -0.001 0.000 

0.026 -0.002 0.003 0.000 0.000 

0.004 -0.001 0.000 0.000 0.000 

0.000 0.000 0.000 0.000 0.000. 

0.576 -0.121 0.048 0.006 0.000 

-0.121 0.164 0.020 -0.002 -0.003 

0.048 0.020 0.010 0.000 -0.001 

0.006 -0.002 0.000 0.000 0.000 

0.000 -0.003 -0.001 0.000 0.000 



Finally, the probabilities are 



Pc = 



0.80703 0.08622 0.02034 0.08622 

0.08622 0.80703 0.08622 0.02034 

0.02034 0.08622 0.80703 0.08622 

0.08622 0.02034 0.08622 0.80703 



: 0.80703 



: 0.19297 . 



B. Performance of 4-PSK and 8-PSK 

The SRM-GUS approach has been applied to evaluate the 
error probability Pg in 4-PSK and 8-PSK systems following 
the steps outlined above. The results are illustrated in Fig. 3 




Error probability in 4-PSK versus Ns for some values of 
K. Solid lines refer to quantum detection and dashed lines 
to classical homodyne detection. 



and in Fig. 4, where Pg is plotted versus the average number of 
photons per symbol Ng for some values of the thermal noise 
parameter >[. In particular, the curve for 3V = 0, which refers 
to the absence of thermal noise, was checked with the results 
obtained with pure states (see (18)) and a perfect agreement 
has been found. To assure an overall accuracy of e = 10~^, 
the reference density operator po was approximated with a 
matrix of size = 145 with a rank running from 1 to 48 in 
dependence of 'N. 

In Fig. 3 and Fig. 4 the PSK quantum detection is compared 
with homodyne counterpart. Remarks similar to that made for 
16-QAM can be repeated. 
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Fig. 4 — Error probability in 8-PSK versus for some values of 
!N. Solid lines refer to quantum detection and dashed lines 
to classical homodyne detection. 



C. SRM vs Optimal Measurement 

For mixed states, the SRM approach is not optimal, at least 
in general, so that a comparison with the optimal performances 
is adequate. As mentioned above, the maximum probability of 
correct detection is the maximum of X)"=o^ "^^{nPi^i) under 
the constraints that the I!,- are PSD and resolve the identity 
I-K or, equivalently [11], the minimum of Tr(F) under the 
constraint that the operators Y — qiPi are PSD. This is a 
problem of convex semidefinite programming. The numerical 
evaluation of the optimum for the applications considered 
in the paper have been performed by Matlab with the LMI 
(Linear Matrix Inequality) Toolbox. A comparison is presented 
in Fig. 5 for the 2-PSK and 4-PSK systems and shows that 
the SRM exhibits an error probability about 30% greater than 



10-* 



10"^ 









2-PSK 








SRIN 


[ 














N=0.2 




( 

(I 


)pumum N> 
lelstrom) 

















1 2 3 4 5 

N, 

Fig. 5 - Error probability in 2-PSK and 4-PSK versus N, 
the optimum detection. 



the optimum detector. Since similar results hold true also for 
8-PSK and 16-QAM, we conclude that SRM is "pretty good" 
also in the presence of thermal noise. 

Remark. For the 2-PSK (as for any other binary format) the 
exact evaluation can be carried out using Helstrom's theory 
(see [10]). This possibiUty was used to check the results 
obtained with the Matlab LMI toolbox. 

VII. Conclusions 

We calculated the error probabiUty of QAM and PSK quantum 

data communication system in the presence of thermal noise 
with quantum detection based on the SRM technique. The 
main novelty of the paper lies in the performance evaluation of 
quantum data transmission system affected by thermal noise, 
not necessarily in a small amount. The lack of results in the 
literature about this topic is surely due to the difficulties of 
numerical computation of optimal detection. The extension 
of the SRM approach to mixed states by Eldar et al. [6] 
allowed us to develop such computations with a relatively 
limited amount of numerical complexity. 

Comparisons made with the performance of classical ho- 
modyne detection give results similar to that for OOK and 
BPSK schemes and evidence the superiority of the quantum 
detection also in the presence of thermal noise. A comparison 
of the SRM with the optimal detection performance, evaluated 
by a convex semidefinite programming package, shows only 
a moderate impairment, so that the obtained results can be 
considered a very good approximation of the optimal perfor- 
mances. 

The results of the application of the SRM to geometrically 
uniform symmetric states will enable one to consider other 
quantum modulation schemes, both in absence and in presence 
of thermal noise. In particular, pulse position modulation 
(PPM) has recently been considered [3] for possible appU- 
cations to deep space quantum communications. 



10"^ 



10-^ 









4-PSK 










SRIV 


- 

2^^=0.4 








)ptimurn\ 


5sr=o.2 ^ 













for some values of K. The SRM detection is compared with 
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However, further problems of computational complexity 
arise, owing to the fact that the natural model for the PPM 
scheme is the tensorial product of Hilbert spaces. Our research 
on the topic is in progress. 
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